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The general solution of SUSY intertwining relations for three-dimensional 
Schrodinger operators is built using the class of second order supercharges 
with nondegenerate constant metric. This solution includes several mod- 
els with arbitrary parameters. We are interested only in quantum systems 
which are not amenable to separation of variables, i.e. can not be reduced to 
lower dimensional problems. All constructed Hamiltonians are partially in- 
tegrable - each of them commutes with a symmetry operator of fourth order 
in momenta. The same models can be considered also for complex values of 
parameters leading to a class of non-Hermitian isospectral Hamiltonians. 



1 Introduction. 



The method of supersymmetrical (SUSY) intertwining relations appeared naturally in 
the framework of a supersymmetrical approach to nonrelativistic Quantum Mechanics [1]. 
Briefly speaking in the most general formulation the method of SUSY intertwining relations 
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for investigation of different models in Quantum Mechanics means the following [2]. Let us 
try to solve the operator equations^ 



where Hi 2 is a pair of Schrodinger Hamiltonians, and the intertwining operator Q + is called 
the supercharge. This relation and its Hermitian conjugate with supercharge Q~ = (Q + Y 
lead to the isospectrality of Hamiltonians Hi j2 up to possible zero modes of supercharges. 
The bound state eigenfunctions of Hi 2 are related (up to normalization factors) by the 
supercharges: 



¥^(x) = Q~¥^(x); = Q+^)(£) ; x = (x l7 x 2 ,x 3 ); n = 0,1,2,... (3) 



If Q + has some zero modes, and they coincide with the wave functions of H 2 (or corre- 
spondingly, Q~ has some zero modes coinciding with the wave functions of Hi), these wave 
functions are annihilated according to (2) and have no analogous states in the spectra of the 
partner Hamiltonian. 

This scheme works independently on the nature and specific properties of operators Hi 2 . 
In particular, it was realized in one-dimensional space with of first order in derivatives 
both for scalar [1] and matrix [3] potentials in the stationary formulation of Schrodinger equa- 
tion. The case of nonstationary Schrodinger equation was studied in [4]. The SUSY inter- 
twining relations with second order supercharges Q ± were introduced first in one-dimensional 
case in [5], including new irreducible supercharges which cannot be factorized with two first 
order multipliers. 

A very useful role was played by second order supercharges in the case of two-dimensional 

SUSY intertwining relations [6]. A class of solutions of such intertwining relations was 

obtained. Each of these Hamiltonians is completely integrable, commuting with its nontrivial 

d From a pure mathematical point of view, the problem can be formulated as follows. We look for two 
different factorizations of the operator of fourth order in derivatives (1) such that H\ t 2 are of Schrodinger 
form, and Q + is the same intertwining second order multiplier. 



HiQ + = Q + H 2 ; 



(1) 




(2) 
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symmetry operator of fourth order in momenta R\ = Q + Q~ or R 2 = Q~Q + . Sometimes, 
besides being completely integrable the model is partially solvable: part of its spectrum and 
corresponding wave functions are constructed by special methods like shape invariance and 
SUSY separation of variables [7], [6]. In the case of two-dimensional generalization of Morse 
potential the exact solvability of the model was proven [8]. It is necessary to stress that only 
models, which do not allow for conventional separation of variables, were considered in the 
papers on two-dimensional Quantum Mechanics mentioned above. 

Some attempts were made to generalize the method of SUSY intertwining relations to the 
physically most interesting case of three-dimensional space. First, the direct generalization 
of formalism with the first order supercharges was constructed for arbitrary dimensionality 
of space in [2]. This approach has a specific property: the intertwining relations link scalar 
Hamiltonians with a chain of matrix Schrodinger-like operators. These Hamiltonians can be 
interpreted as operators describing quantum systems with some internal degrees of freedom 
(for example, spin [9]), but rather frequently this property seems to be inconvenient. Inter- 
twining of two scalar Hamiltonians by first order supercharges was shown [10] to offer only 
solutions with conventional separation of variables. The only alternative generalization [11], 
to our knowledge, is analogous to the idea elaborated earlier in two-dimensional case [6]: to 
use the intertwining relations between pairs of scalar Hamiltonians but with second order 
supercharges. In paper [11] the particular solutions of such intertwining relations were found 
for the special form of second order supercharges ^f^- — + .... 

In the present paper we will consider the most general form of supercharges Q ± with 
constant metric g^, i.e. Q ± = Qik-^rr^ + where the summation over repeated indices 
i, k — 1, 2, 3 is implied here and below. In contrast to the similar approach in two-dimensional 
space [6], where a list of particular solutions was found, in three-dimensional case the general 
solution of intertwining relations will be built for arbitrary constant metric excluding 
the degenerate (up to rotations) case of g^ = (1,0,0). This difference is due to a more 
restrictive character of intertwing relations in higher dimensions. The general solution of 
three-dimensional intertwining with constant metric in supercharge includes several classes 



3 



of models, which depend on some parameters. These parameters can be chosen either real 
or complex leading to Hermitian or non-Hermitian partner Hamiltonians. 



2 General solution of second order intertwining rela- 
tions with nondegenerate metric. 

As announced in Introduction, we shall study the intertwining (1) with the most general 
supercharges of second order in derivatives with constant metric (highest order coefficients) 

Q + = g l Ad k + C l {x)d l + B{x)- d t = ^. (4) 

By space rotations the matrix g ik can be diagonalized to = gaS^. At first, we will be 
interested in situations when all three diagonal elements gu after such rotation do not vanish. 
The case gu = (1,-1,0) was considered in [11], the case gu = (1,1,0) will be considered 
later in this Section, the case gu = (l,a, 0); a ^ ±1 gives no interesting solutions for (1), 
and for the degenerate case gu = (1, 0, 0) only some particular solutions will be given in the 
very end of this Section. 

Thus, let us consider now the supercharges Q ± , which by a suitable normalization can 



be reduced to two classes of gi k : 

(A) 9ii = (l,l,d); d^Od^l (5) 

(B) gu = (l,d 2 ,d 3 y, 4^0,1; d 2 ^ 0,1; d 2 ^ d 3 . (6) 

The intertwining relations (1) for constant g ik can be rewritten as a system of 6+3+1 = 10 
differential equations, some of them nonlinear: 

d l C k (x) + d k C l (x) = 2V(x)g tk ; (7) 

A^d(x) + 2d l B(x) + 2g ik d k V 2 = 2V(x)C l (x); (8) 

A^B(x) + g lk O l d k V 2 {x) + C t (x)d z V2(x) = 2V(x)B(x); (9) 

Vx{x) - V 2 {x) = 2V(x). (10) 
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The six equations (7) and the defining Eq.(lO) can be combined: 

&Ci(£) = d^d 2 C 2 (x) = d^d 3 C 3 (x) = V(x); (11) 
d z C k (x) + d k C l (x)=0.; i^k. (12) 

After further differentiation of (12) in respect to Xj with j ^ i, j ^ k we obtain: 

Ci(x) = U{x u x 2 ) + f 2 (xi,x 3 ); 
C 2 (x) = [gi(x 2 ,x 3 ) + g 3 (x 1 ,x 2 )]d 2 ; 
C 3 (x) = [h 1 (x 2 ,x 3 ) + h 2 (x 1 ,x 3 )]d 3 , 

where f,g,h are auxiliary functions specified in the following. 

By substitution of these relations into (11), we get the functional-differential equations, 
which are integrated in a general form: 

gi(x 2 ,x 3 ) = x 2 G' 3 (x 3 ) + G 2 (x 2 ) + K 3 (x 3 )\ 
f 2 (x 2 ,x 3 ) = x 1 G' 3 {x 3 ) + G 1 (x 1 ) + M 3 {x 3 ); 
f 3 (x 1 ,x 2 ) = Gi(xi) + xiG' 2 (x 2 ) + K' 2 (x 2 ) - Gi(xi); 
g 3 (x 1 ,x 2 ) = G 2 (x 2 ) + x 2 G' 1 (x 2 ) + K[(xi) - G 2 (x 2 ); 
h 2 (xi,x 3 ) = x^'^xi) + Li{xi) + M 3 (x 3 ); 
h 1 (x 2 ,x 3 ) = x 3 G' 2 (x 2 ) + L 2 (x 2 ) + G 3 (x 3 ) - M 3 (x 3 ). 

Therefore, we have the expressions for coefficient functions Ci{x) of the supercharges in 
terms of the eleven arbitrary functions above Gi(xi), Gi(xi), G 2 (x 2 ).... If these expressions are 
inserted into equations (12), the eleven functions above will be specified being polynomials 
maximally of second order in space coordinates. Thus the general solution of (11) takes the 
form: 

Ci(x) = ai(x\- d 2 x\ - d 3 xl) + 2a 2 xix 2 + 2a 3 x\x 3 + bxi -bix 2 -b 3 x 3 + Ki\ 

C 2 (x) = a 2 (-xj + d 2 x\ - d 3 x 2 3 ) + 2a 1 d 2 x 1 x 2 + 2a 3 d 2 x 2 x 3 + bd 2 x 2 + bxXi - b 2 x 3 + k03) 

C 3 (x) = a 3 (-x\ - d 2 x\ + ^3X3) + 2a 1 d 3 x 1 x 3 + 2a 2 d 3 x 2 x 3 + ^3^3 + hxi + b 2 x 2 + k 3 \ 



and the difference between potentials Vi,2 is given now by: 



V(x) = 2(a 1 x 1 + a 2 x 2 + a 3 x 3 ) + b. 



(14) 



The next step consists in solving the three equations (8), where the first Laplacian terms 
are constants defined from (13): 



Now we will restrict ourselves to the case of metric (5), i.e. d 2 — 1, d 3 = d ^ 0, 1. Then, 
after the simple manipulations (derivatives and linear combinations) with (15) we obtain the 
necessary conditions: 



As for the constant a 3 , it is convenient to consider separately two options: (i) a 3 ^ and 
(ii) a 3 = 0. 

In the case (i) a 3 / 0a suitable translation of space coordinates allows to cancel the 
linear terms in (13): 

Ci(x) = 2a 3 x 1 x 3 + ei, C 2 (x) = 2a 3 x 2 x 3 + e 2 ; C 3 (x) = a 3 (dx 2 3 -x\- x\) + e 3 ; V(x) = 2a 3 x 3 , 



where e« are new arbitrary constants. Taking into account the numerical values of the 
constants, Eq.(16), we obtain for the considered case (d 2 = 1; a 3 ^ 0) the general formulas: 



diB(x) + diV 2 (x) = V(x)Ci(x) - ci; ci = a ± (l - d 2 - d 3 ); 
d 2 B(x) + d 2 d 2 V 2 (x) = V(x)C 2 (x) - c 2 ; c 2 = a 2 (d 2 - 1 - d 3 ); 
d 3 B(x) + d 3 d 3 V 2 (x) = V(x)C 3 (x) - c 3 ; c 3 = a 3 (d 3 - 1 - d 2 ). 



(15) 



cl\ = a 2 = b\ = 0. 



(16) 



(17) 



V 2 (x) 



U, 3 U 4 x ro 2 2/ 2 , 2\ 2 , , 

= 2 (d- l) ^ 3 ~ ^ITl^Wi + z 2 ) - a 3 e 3 a; 3 + c 3 a; 3 + 
+ 2a 3 x 3 (eiXi + e 2 x 2 ) + g 3 (a; 3 ) - F 3 (x 1 ,x 2 )]; 




(18) 



= ~ 2 (^_ + rf3i"^ 2rf + 1 ) a 3 a; 3( a; i + ^2) - a 3 e 3 a; 3 + c 3 x 3 + 
+ 2rfa 3 a; 3 (eiXi + e 2 x 2 ) + +<ig 3 (a; 3 ) - F 3 (xi,x 2 )], 




(19) 
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where q 3 (x 3 ) and F 3 (xi,x 2 ) are arbitrary functions, which will be defined from the last 
intertwining relations Eq.(9). 

We use the divergence of vector equation (8) to transform (9) to: 

-(1 + ^)A ( VV(x) + Ci{g)diV(g) + (2 + d)V 2 (x) + Q(x)d t V2(x) - 2V(x)B(x) = 0. (20) 

Since the l.h.s. can be written as a polynomial in x 3 , we put the coefficients to zero. Insert- 
ing explicit expressions (17), (18), (19) into (20) and after rather long but straightforward 
algebra, we derive that the solution exists only for the metric (1, 1, —1), i.e. for d = —1. In 
this case the function q 3 is linear (with arbitrary constant «o): 

^3(^3) = -CL3X3 + «o- (21) 

As for the function F 3 , it is defined by two differential equations: 

-2(F 3 (xi, x 2 ) + oq) - (andi + x 2 d 2 )F 3 (x 1 ,x 2 ) + (e 2 + e\) + (e 3 - a 3 p 2 )(3a 3 p 2 - e 3 ) = 0; (22) 

(eidi + e 2 d 2 )F 3 (x 1 ,x 2 ) = 2a 3 (e 3 - a 3 p 2 )(eiXi + e 2 x 2 ), (23) 

where we introduced the radial coordinate in the plane p 2 = x\ + x\. 

In the particular case of constants t\ = e 2 = the equation (23) is trivially satisfied, and 
the general solution for F 3 can be represented in polar coordinates p, as: 

F 3 ( Xl ,x 2 ) = -\a\p* + a 3 e 3 p 2 - - I(2« + e 2 3 ). (24) 

Thus we obtain the expressions for the partner potentials: 

Vl , 2 (x) = £(xl + p 4 ) + ^fxtp 2 - a -^{xl + p 2 ) ± 2a 3 x 3 + M + Const, (25) 

which are anharmonic oscillators of fourth order with additional 1/p 2 term. 

In the particular case e x e 2 7^ it is convenient to introduce the new space coordinates: 



so that in (22), (23) 



ei<9i + e 2 d 2 = 2d yi ; + x 2 d 2 = yid yi + y 2 d y2 . (27) 

Then the equation (23) can be solved in a general form: 

aHH + 3el - Hel) y2y2 + + ^ + 2(e? _ + (2g) 



16 1 ^ 4 

with a function p 2 (y 2 ) defined from (22): 



a 3 e 3 (e? + ej) 2 e| - ef - e| + 2q « o|(e| + e 2 .) 2 4 

p 2 (y 2 ) i y 2 - 2 + ^ 32 y 2 . 



Now, inserting this function into (28) and going back to the variables Xi, we obtain the final 
expression for potentials Vi >2 : 

Vi, 2 {x) = ^(x A 3 + p 4 ) + ^4° 2 - ^f%3 + P 2 ) + «3e3(eiX! + e 2 x 2 ) - 

2 2 

M ? ] ± 2a 3 x 3 + Const. (30) 



2(e 2 xi - e\_x 2 ) 2 

The last option for the values of e±, e 2 we have to consider is the case with one of them 
vanishing and the other not. For definiteness, let us take e 2 = 0. Then the general solution 
of (23) is: 

F 3 (x l7 x 2 ) = yx 2 (2e 3 - a 3 x\ - 2a 3 x 2 2 ) + q 2 (x 2 ), (31) 

with an arbitrary function q 2 (x 2 ). Using it in (22), we obtain a first order differential equation 
for q 2 (x 2 ), which admits the general solution: 

Q2{x 2 ) = -\a\x\ + a 3 e 3 x\ + -j + \{e\ - e\ - 2a ) (32) 
z x 2 z 

(p is an arbitrary integration constant). The resulting potential has the form: 

V h2 (x) = ^{x\ + p 4 ) + ^x 2 3 p 2 -^(x 2 3 + p 2 ) + a 3 e lXl x 3 ± 2a 3 x 3 - 4 + Const (33) 
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In the case (ii) a 3 = equations (13), (16) together with a suitable shift of space coordi- 
nates allow to conclude that the coefficient functions Cj(x) become linear functions: 

C\(x) = bxi - b 3 x 3 ; C 2 (x) = bx 2 - b 2 x 3 ; C 3 (x) = bdx 3 + b 3 x 1 + b 2 x 2 ; V(x) = b (34) 

with arbitrary real constants b,b 2 ,b 3 . Then from the explicit formulas for Vi,2 we find that 
both partner Hamiltonians Hi t2 being nonisotropic harmonic oscillator are amenable to sep- 
aration of variables. 

But there is just one solution for the case a 3 = which does not allow separation 
of variables. In order to obtain it, we have to extend the original assumption about the 
constants in (34). Indeed, if we allow b 2 and b 3 to be complex numbers with b 2 = ib 3 , after 
simple calculations we obtain: 



with arbitrary constants n, u. 

Both potentials are second order polynomials with complex coefficients 6 and they differ 
from each other by the arbitrary constant b. Such pair of intertwined Hamiltonians realize the 
simplest kind of shape invariance [14], [7]. In this context, it means that if the Hamiltonian 



an arbitrary eigenvalue Eq in the spectra of Hi is accompanied by the tower of eigenstates 
with energies E + 2bk. The same is obviously true for H 2 . We must stress that this property 
of the spectra are formulated up to possible zero modes of supercharges Q ± . These zero 
modes can, in principle, truncate the tower described above. 

An additional property of the spectrum of model (35) is provided by its symmetry. In- 
deed, for real values of parameter n both Hamiltonians H\ 2 are invariant under the combined 

e The complex coefficients appear in the supercharges Q also. See detail discussion of one-dimensional 
case in [12] and of two-dimensional case in [13]. 




(35) 




(36) 
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reflection P 2 , i.e. x 2 — > —^2, and time reflection T. Due to this P 2 r-symmetry, the spectrum 
of Hamiltonians consists of [15], [16] real eigenvalues E = E* and complex conjugate pairs 
E, E*. 

One can notice that all models above which were supposed to have real potentials can 
be easily complexified also by chosing complex values of parameters in (25), (30), (33). 
Although no property of shape invariance is observed for these potentials, some invariances 
under combined reflections can be realized with suitable choice of complex parameters. For 
example, P{T— symmetry of potential (30) for pure imaginary e x and other parameters real. 

Let us consider now the case when the metric can be reduced to 

(C) 9ik = (1,1,0). (37) 

Then the general solution of Eq.(7) is: 

Ci = -x 3 (2axi + an) - 2{vx\ + u ± ); C 2 = -x 3 (2ax 2 + a 2 ) - 2{vx 2 + u 2 ); 

C 3 = a(x\ + x\) + aiXi + a 2 x 2 + (3; V = -2(ax 3 + u), (38) 

with arbitrary constants a, a lj2 , u, v 12 . For a 7^ 0, it can be transformed by translations to: 

Ci = -2(axxx 3 + z/i), Ci = -2(ax 2 x 3 + v 2 )\ C 3 = a(x\ + x\) + (3; V = -2ax 3 . 

The first two equations of Eq. (8) lead to: 

B + V 2 — 2ax 3 [a(xl + xl)x 3 + 2{y x x x + v 2 x 2 )\ + f 3 (x 3 ), (39) 

with an arbitrary function f 3 (x 3 ). The last equation of Eq.(8) gives B(x): 

B(x) = -axl[a(x\ + x\) +13}- 2ax 3 + F 3 (x 1 ,x 2 ), (40) 

where F 3 (xi,x 2 ) is also an arbitrary function, and (3 is a constant. Therefore the potential 
V 2 (x) is obtained: 

V 2 (x) = 3a 2 xjp 2 + apxj + 2a{2v lXl + 2v 2 x 2 + l)x 3 + f 3 - F 3 . 
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The arbitrary functions f 3 (x 3 ), F 3 (xi,x 2 ) are defined by Eq.(9) using also Eq.(8): 

f 3 (x 3 ) = 4a 2 x 4 3 + fox 3 3 + I3 2 xl + I3 x x 3 + fo; 

j? , \ 1 2 4 1 n 2 f(<f>) 

F 3 {x 1 ,x 2 ) = --a p - -fop a~, 

2 2 p z 

with arbitrary periodical function f((f>) of polar coordinate 0. The constants above must be 
chosen vanishing: v\ — v 2 — fo — fo — 0. Finally, the partner potentials V 1}2 {x) are obtained 
as: 



V 1:2 = a 2 (3x 2 3 p 2 + 4xt + -p 4 ) + fo{xi + -p 2 ) + 2ax 3 . (41) 

In the case a — 0, Eq.(38) gives after translation: 

Ci = -(0:1X3 + IviXi); C 2 = -(a 2 x 3 + 2ux 2 ); C 3 = oi^i + o 2 x 2 ; V = -2u, (42) 
and, using (8): 

B = -2u(a±xi + a 2 x 2 )x 3 + A^ 3 (xi, x 2 ); (43) 

V 2 = 2v 2 p 2 + 4i/(oiXi + 02^2) + n 3 (x 3 ) - N 3 (x u x 2 ), (44) 

with arbitrary functions n 3 (x 3 ) : N 3 (x± : x 2 ). The analysis of Eqs.(42)-(44) together with 
Eqs.(8), (9) shows that the only solution with real parameters corresponds to models al- 
lowing standard separation of variables, which are not studied in this paper. 

The only constant metric g^ in Q ± , which does not allow to find the general solution, 
is (up to rotations) the degenerate one - gu = (1, 0, 0). In this case only particular solutions 
were found: 

t r 1/ . \2 -2/3 5C 2 2/3 . 16C 2 2 

V 2 = -{ax 2 + x 3 ) x ± - — {ax 2 + x 3 ) x{ +7^7— — ¥ r{ax 2 + x 3 ) - 
9 6 9(1 + cr) 



a 4 c z 2 I + a 2 4/3 3a 4 2 /3 7 



ax 2 + x 3 ) + — x 1 H — Xi + —Xi + — x x +&(ax 3 -x 2 ), 



l + o 2V 4 1 2 1 8c 1 36 

V = -(ax 2 + x 3 )xi 4/3 + c 
9 
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and 



V 2 = x\[?>(ax 2 + x 3 ) 2 - 3b(ax 2 + x 3 ) - 2(b 2 + c)] + 



1 + a 2 




x\ — UJX 1 2 + 



1 



-[A(ax 2 + x 3 ) 4 - 8b(ax 2 + x 3 ) 3 - c(ax 2 + x 3 ) 2 + 



1 + a 



+b(4b 2 + c)(ax 2 + ar 3 )] + <&(ax 3 - x 2 ), 
V = -2(ax 2 + x 3 ) + b, 

where $ is an arbitrary function of its argument. 

3 Conclusions. 

The aim of the present paper was to provide an exhaustive analysis of intertwining re- 
lations for second order supercharges with constant metric in three-dimensional case. In 
this space the intertwining relations are rather restrictive: only specific classes of solutions 
exist. Namely, four different classes of real potentials were built in the framework of general 
solution of intertwining relations with nondegenerate metric - (25), (30), (33), (41). All of 
them have the form of anharmonic oscillator of fourth order in coordinates plus the term 
with 1/x 2 — dependence. In addition, particular solutions were obtained for degenerate metric 
(1,0,0). For all these models pairs of almost (up to zero modes of isospectral Hamil- 
tonians were obtained. Each model is at least partially integrable, since the Hamiltonians 
commute with the corresponding symmetry operators of fourth order in momenta. 

Although the intertwining relations themselves do not provide, in general, the wave func- 
tions and energy eigenvalues, in two-dimensional space the specific SUSY methods - SUSY 
separation of variables and shape invariance [7], [6], [8] - allowed to find a part or even 
the whole spectrum. In contrast to this situation, it seems to be impossible to apply these 
methods in three-dimensional space: there is no shape invariance property and there are no 
tools to find zero modes of Q ± explicitly. 

It was also useful to consider the obtained results from the point of view of Quantum 
Mechanics with non-Hermitian Hamiltonians [12], [13], [15], [16]. In this case, the potentials 
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with complex values of parameters obey some discrete symmetries with antilinear operators, 
leading to corresponding specific properties of complex spectra. 
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